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the three-dimensional coordinate- mesh representation, we suggest that, 
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1 Introduction 



Recently, super deformed (SD) rotational bands have been discovered in 36 Ar, 
40 Ca and 44 Ti jUEHSHHEHEl- One of the interesting new features of them 
is that they are built on excited + states and observed up to high spin, in 
contrast to the SD bands in heavier mass regions where low-spin portions of 
them are unknown in almost all cases [301311011111 • These excited + states 
may be associated with multiparticle-multihole excitations from the spherical 
closed shells, so that we can hope to learn from such data detailed relationships 
between spherical shell model and SD configurations. For the mass A =30-50 
region, although existence of a SD band in 32 S with the SD magic number 
N = Z = 1Q has been expected for a long time J2], it has not yet been 
observed and remains as a great challenge [T3 l lT H IT5 | ll6 [ IT7j . 

In this paper, as a continuation of the systematic theoretical search P~^.fTH| 
for SD bands in the mass A =30-50 region by means of the cranked Skyrme- 
Hartree-Fock (SHF) method jTH], we would like to suggest that, in addition 
to 32 S, the neutron-rich nucleus 36 S and the nuclei, 48 S and 50 S, which are 
situated close to the neutron-drip line [2D1I2I], are also good candidates for 
finding SD rotational bands in Sulfur isotopes. The appearance of the SD 
band in 36 S is suggested in connection with the SD shell structure at N = 20 
characterizing the observed SD band in 40 Ca. The drip-line nuclei, 48 S and 
50 S, are expected to constitute a new "SD doubly closed" region associated 
with the SD magic numbers, Z = 16 for protons and iV ~ 32 for neutrons. An 
interesting theoretical subject for the SD bands in nuclei near the neutron drip 
line is to understand deformation properties of neutron skins. The calculated 
density distributions indeed exhibit superdeformed neutron skins. 

The calculation has been carried out with the use of the three-dimensional 
(3D), Cartesian coordinate- mesh representation without imposing any sym- 
metry restriction ^31E3- in parallel, we also carry out the standard calcula- 
tions [221 12H1I2H 1213 imposing reflection symmetries. By comparing the sym- 
metry restricted and unrestricted calculations, we can examine the stability of 
the SD solutions of the SHF equations against reflection-asymmetric deforma- 
tions. In this way, we have found several cases where the SD minima obtained 
in the symmetry-restricted calculations are in fact unstable with respect to 
the reflection-asymmetric deformations. In general, the SD states are rather 
soft against reflection-asymmetric deformations, so that we need careful study 
about the stabilities of them against various kinds of deformation breaking the 
reflection symmetries. 

After a brief account of the cranked SHF calculational procedure in Sec- 
tion 2, we present and discuss results of the calculation in Section 3, and give 
conclusions in Section 4. We shall present deformation energy curves for Sulfur 
isotopes from 32 S to 50 S, and focus our attention on properties of rotational 
bands built on the SD + states, stabilities of the SD local minima against 
the reflection-asymmetric deformations, and density distributions of the SD 
states. 
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A preliminary version of this work was reported in [*""""[ . 

2 Cranked SHF calculation 

Since the cranked SHF method in the 3D coordinate-mesh representation is 
well known fl^\ 123*1 "2*"""( l25j. we here give only a minimum description about 
the computational procedure actually adopted. For a recent comprehensive 
review on selfconsistent mean-field models for nuclear structure, see Ref. [191 . 
The cranked SHF equation for a system uniformly rotating about the X~£lXlS 
is given by 

5 <H - u rot J x >= 0, (1) 

where H, u rot and J x mean the Hamiltonian with the Skyrme interaction, the 
rotational frequency and the x-component of angular momentum, respectively, 
and the bracket denotes the expectation value with respect to a Slater determi- 
nantal state. We solve the cranked SHF equation by means of the imaginary- 
time evolution technique [22] in the 3D Cartesian-mesh representation. The 
algorithm of numerical calculation is the standard one [2"2^l2ll[l2""""( l25|. except 
that we allow for both reflection- and axial-symmetry breakings. In this case, 
it is important to accurately fulfill the center-of-mass and principal-axis condi- 
tions. This is done by means of the constrained HF procedure [2Zj . We solved 
these equations inside the sphere with radius R=10 fm and mesh size h=l fm, 
starting with various initial configurations. The accuracy for evaluating de- 
formation energies with this mesh size was carefully checked by Tajima [23] 
and was found to be quite satisfactory. When we make a detailed analysis of 
density distributions, however, we use a smaller mesh size of h = | fm. In 
addition to the symmetry-unrestricted cranked SHF calculation, we also carry 
out symmetry-restricted calculations imposing reflection symmetries about the 
(x, y)-, (y, z)- and (z, a;)-planes. Below we call these symmetry-unrestricted 
and -restricted cranked SHF versions "unrestricted" and "restricted" ones, re- 
spectively. 

Solutions of the cranked SHF equation give local minima in the deforma- 
tion energy surface. In order to explore the deformation energy surface around 
these minima and draw deformation energy curves as functions of deformation 
parameters, we carry out the constrained HF procedure with relevant con- 
straining operators [2Z] • For the Skyrme interaction, we adopt the widely used 
three versions; SIII [28 , SkM* PSJ and SLy4 fSH]. 

3 Results of calculation 

3.1 Deformation energy curves 

Figure 1 shows deformation energy curves for Sulfur isotopes from 32 S to 50 S 
obtained with the use of the SIII interaction. Solid lines with and without 
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filled circles represent the results obtained by the unrestricted and restricted 
versions, respectively. The result of calculation indicates that the SD minima 
(with the quadrupole deformation parameter fa ~ 0.6) appear in the neutron- 
rich nucleus 36 S and the drip-line nuclei, 48 S and 50 S, in addition to the well- 
known case of 32 S. As seen in Figs. 2 and 3, similar results are obtained for 
the SkM* and SLy4 interactions, except that the SD states in 48 S is unstable 
against the reflection-asymmetric deformation for the SLy4 interaction (see 
subsection 3.3). 

As discussed in Refs. |T4* | ll5 p i6 [ lT7 t ll8| . the SD local minimum in 32 S corre- 
sponds to the doubly closed shell configuration with respect to the SD magic 
number Z = N = 16 and involves two protons and two neutrons in the down- 
sloping single-particle levels originating from the fa/2 shell. The SD local 
minimum in 36 S results from the coherent combination of the SD magic num- 
ber, Z = 16, and the neutron shell effects occurring at large deformation for 
N = 20. The latter shell effect has been confirmed recently by the discovery of 
the SD rotation band in 40 Ca |U|3]. The SD shell gap at iV = 20 is associated 
with the 4p-4h excitation from below the N = 20 spherical closed shell to the 
fa/2 shell. Focusing our attention on the occupation numbers of such high-j 
shells and distinguishing protons(7r) and neutrons(zz), these SD configurations 
in 32 S and 36 S are denoted in Figs. 1-3 as J 2 / 2 and J 2 / 4 , respectively. 

The SD local minima in the drip- line nuclei, 48 S and 50 S, result from the 
coherent combination of the proton SD shell effect and the neutron shell effects 
occurring at superdeformation for N =32-34. The neutron configurations in 
these SD states are similar to those in the known SD bands in 60 Zn and 62 Zn 
associated with the SD magic numbers N =30-32 |3T|l32j . We find that the SD 
shell effect is strong also for N = 34 in the Sulfur isotopes under consideration, 
while the SD local minimum in 46 S with N = 30 is unstable against the 
reflection-asymmetric deformation (see subsection 3.3). In the drip-line nuclei 
48 S and 50 S, the fa/2 shell is fully occupied even in the spherical limit and the SD 
configurations involve neutron excitations from the /p-shell to the g$/2 shell. 
As before, focusing our attention on the occupation numbers of the high-j 
shells, let us use the notation f^g™ 2 for a configuration in which single-particle 
levels originating from the fa/ 2 and gg/ 2 shell are occupied by n\ protons and 
n,2 neutrons, respectively. With such notations, the SD local minima in 48 S 
and 50 S correspond to the / 2 g 2 , and / 2 g 4 configurations, respectively. 

The appearance of the SD minimum in 36 S suggests that we can expect a SD 
band associated with the same neutron configuration to appear also in the iV = 
20 isotone, 38 Ar, situated between 36 S and 40 Ca. We examined this point and 
the result is shown in Fig. 4. We find that the two local minima associated with 
the configurations / 2 / 2 and / 2 / 4 compete in energy and their relative energy 
differs for different versions of the Skyrme interaction: As clearly seen in the 
deformation-energy curves obtained by the symmetry-restricted calculations, 
the former with smaller fa is slightly lower for SkM* and SLy4 while the 
latter with larger fa is slightly lower for SIII. Counting both protons and 
neutrons, these local minima respectively correspond to the 4p-6h and 6p-8h 
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configurations with respect to the spherical doubly closed shell of 40 Ca. As 
we discussed in the previous papers [HH^H], the two configurations can mix 
each other in the crossing region through the reflection-symmetry breaking 
components in the mean field. Specifically, around the crossing point between 
the down-sloping [321 1] level (coming from the shell) and the up-sloping 
[200 i] level (coming from the shell below the N = 20 spherical magic 
number), the r 3 Y"3i-type non-axial octupole deformation is generated, and they 
mix each other through this component of the mean field (see Fig. 5). Note 
that the matrix element of the operator r 3 Y" 31 between the two levels satisfies 
the selection rules, An 3 = 2 and A A = 1, for the asymptotic quantum numbers 
713 and A. As a result of this mixing, the deformation-energy curve becomes 
rather flat in the symmetry-unrestricted calculation. Recently, the SD band 
corresponding to the 4p-6h configuration was found in experiment [33]. The 
data suggest significant competition between different configurations, which 
requires further analysis of shape fluctuation dynamics by going beyond the 
static mean-field approximation. 

3.2 SD rotational bands 

Let us focus our attention on the SD local minima shown in Figs. 1-3, and 
investigate properties of the rotational bands built on them. Excitation ener- 
gies of these SD rotational bands are plotted in Fig. 6 as functions of angular 
momentum. These rotational bands are calculated by cranking each SHF so- 
lution (the SD local minima in Figs. 1-3) and following the same configuration 
with increasing value of u rot until the point where we cannot clearly identify 
the continuation of the same configuration any more. Thus, the highest values 
of angular momentum in this figure does not necessarily indicate the band- 
termination points but merely suggest that drastic changes in their microscopic 
structure take place around there. Different slopes with respect to the angular 
momentum between 36 S and 50 S can be easily understood in terms of the well 
known scaling factor A 5 ^ 3 for the rigid-body moment of inertia. This point 
can be confirmed in Fig. 7 which displays the angular momentum J, the kine- 
matical and dynamical moments of inertia, = I /uj TOt and = dl /duj TO t, 
and the rigid-body moments of inertia J„ g — m j p(r)(y 2 + z 2 )dr as functions 
of the rotational frequency c<j rot . We see that the calculated moments of in- 
ertia are slightly larger than the rigid-body values at co> rot = 0, and smoothly 
decrease as u rot increases until u TOt «2.5 and 1.8 MeV/h for 32 > 36 S and 50 S, 
respectively. The result calculated with the SLy4 interaction is shown here, 
but we obtained similar results also with the SIII and SkM* interactions. 

Calculated quadrupole deformation parameters (/3 2 ,7) are displayed in the 
upper portion of Fig. 6. We see that the /3 2 values slowly decrease while 
the axial-asymmetry parameters 7 gradually increase with increasing angular 
momentum for all cases of 32 S, 36 S and 50 S. The variations are rather mild in 
the range of angular momentum shown in this figure. Single-particle energy 
diagrams (Routhians) for these SD bands are displayed in Fig. 8 as functions 
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of the rotational frequency uo mt . This figure indicates that level crossings take 
place in 36 S and 50 S if we further increase the angular momentum. 

3.3 Stabilities of the SD states against reflection- asymmetric 
deformations 

Let us examine stabilities of the SD local minimum against both the axially 
symmetric and asymmetric octupole deformation (I30, 131, 132, ^33)- Figure 9 
presents deformation energy curves as functions of the octupole deformation 
parameters (5^ m {m = 0,1,2,3) for fixed quadrupole deformation parameters 
(the equilibrium value of at the SD minimum in each nucleus and 7 = 0). 
The computation was carried out by means of the constrained HF procedure 
with the use of the SIII, SkM*, and SLy4 interactions. The result of calculation 
clearly indicates that the SD states in 32 S, 36 S and 50 S are stable against the 
octupole deformations and that they are softer for j3 3m with lower values of m 
(i.e., for /3 30 and /?3i), irrespective of the Skyrme interactions used. We obtaind 
a similar result also for 48 S (but omitted in this figure). 

Although the SD minima in 32 S, 36 S and 50 S are stable with respect to 
@3m{ m — 0,1,2,3), we found several cases where the SD minima obtained 
in the symmetry-restricted calculations become unstable when we allow for 
reflection-asymmetric deformations of a more general type. As a first example, 
let us discuss the SD minimum in 46 S which appears in the restricted calcula- 
tion (see Figs. 1-3). In this case, the coupling between the down-sloping [330|] 
level (associated with the ft/2 shell) and the up-sloping [202|] level (stemming 
from the d^/2 shell) takes place in the proton configuration, when we allow for 
the breaking of both the axial and reflection symmetries. Thus, the SD config- 
uration flgl mixes with the gl configuration (which lacks the proton excitation 
to the f 7 /2 shell and has a smaller equilibrium value of fo)- As a consequence 
of this mixing, the barrier between the two configurations disappears and the 
SD minimum becomes unstable in the unrestricted calculations (see Figs. 1-3). 
Note that the difference An 3 in the asymptotic quantum number n 3 between 
the two single-particle levels, [330|] and [202|], is three, so that they cannot 
be mixed by the octupole operator r 3 Y" 3 2 which transfers the asymptotic quan- 
tum numbers n 3 and A by An 3 = 1 and AA = 2. Thus, this mixing may be 
associated with the reflection-asymmetric deformation of a more higher order 
like r 5 F 52 . 

As a second example, we take up the SD minimum in 48 S. In this case, 
two configurations, f%gl and f%g^ compete in energy and their relative energy 
differs for different versions of the Skyrme interaction (see Figs. 1-3). When 
we allow for the breaking of both the axial and reflection symmetries, the 
coupling between the down-sloping [43l|] level (associated with the g 9 / 2 shell) 
and the [310^] level in the fp shell takes place in the neutron configuration, 
so that they mix each other. Note that the [431 1] and [310^] levels satisfy 
the selection rules, An 3 = 2 and AA = 1, for the matrix elements of the 
octupole operator r 3 F 31 . In the calculation with the SLy4 interaction, since 
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the former configuration with smaller is situated slightly lower in energy 
than the latter, the barrier between the two configurations disappears as a 
result of this mixing. This mixing effect in conjunction with that mentioned 
above for the f 2 gl configuration in S deteriorates the SD minimum for the 
SLy4 case. 

The above examples indicate detailed microscopic mechanisms within the 
mean-field theory how the stability of the SD local minimum is determined 
by relative energies between the neighboring configurations and their mixing 
properties. 

3.4 Density distributions 

Figure 10 displays the neutron and proton density profiles for the SD states in 
32 S, 36 S and 50 S calculated with the use of the SLy4 interaction. We obtained 
similar results also for SIII and SkM*. In this figure, equi-density lines with 
50% and 1% of the central density in the (x, y)- and (y, z)-planes are drawn 
for the SD bands at / = and at high spins. We can clearly see that superde- 
formed neutron skin appears in 50 S which is situated close to the neutron drip 
line. The root-mean-square values, J (x 2 ), J (y 2 ), yj (z 2 ) and -R rms = J (r 2 ), 
of these density distributions are listed in Table 1. To indicate the deforma- 
tion properties of the neutron skin in 50 S, calculated values for protons and 
neutrons are separately listed together with their sums and differences. We 
obtained density distributions similar to those for 50 S also for the SD state in 
48 S. A similar result of theoretical calculation exhibiting the superdeformed 
neutron skin was previously reported in Ref. |33] for the SD state in the very 
neutron-rich nucleus 2 6gDy 142 . 

4 Conclusions 

On the basis of the cranked SHF calculations in the 3D coordinate-mesh rep- 
resentation, we have suggested that, in addition to the well-known candidate 
32 S, the neutron-rich 36 S and the the drip-line nuclei, 48 S and 50 S, are also good 
candidates for finding SD rotational bands in Sulfur isotopes. Calculated den- 
sity distributions for the SD states in 48 S and 50 S, which are situated close to 
the neutron-drip line, exhibit superdeformed neutron skins. 
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Table 1 



Root-mean-square values, y (x 2 ) , y (y 2 ) , y (z 2 ) and R rms = y(r 2 ), of the 
density distributions at / = (second column) and at / = 20, 22, 28 (third col- 
umn) of the SD band in 32 S, 36 S and 50 S, calculated with the use of the SLy4 
interaction. Neutron and proton contributions are separately listed together 
with their sums (total) and differences (diff.). 
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Figure 1: Deformation energy curves for Sulfur isotopes from 32 S to 50 S calculated 
at / = as functions of the quadrupole deformation by means of the constrained 
SHF procedure with the SIII interaction. The deformation parameter is defined as 
@2 = (X^=i r fY2o(0i, <fii)) / (J2t=i r \ )• The axial-asymmetry parameter 7 is con- 
strained to be zero. Solid curves with and without filled circles represent the results 
obtained by the unrestricted and restricted versions, respectively. The notation 
Z™ 1 /™ 2 indicates a configuration in which single-particle levels originating from the 
f 7 / 2 shell are occupied by ri\ protons and 722 neutrons. Likewise, f™ l g™ 2 indicates 
that levels from the / 7 / 2 shell are occupied by n\ protons and those from the 59/2 
shell by 712 neutrons. 
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Figure 3: The same as Fig.l but for the SLy4 interaction. 
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Figure 4: The same as Fig.l but for 38 Ar and for the SIII, SkM* and SLy4 interac- 
tions. 




Figure 5: Schematic illustration of configuration-mixing mechanism through the oc- 
tupole components of the mean field. When the reflection symmetry is imposed, the 
positive- and negative-parity single-particle levels sharply cross, and the two config- 
urations (having different number of particles in the fj/2 shell) do not mix within 
the mean-field approximation (left-hand side). In contrast, when such symmetry 
restriction is removed, smooth crossover between the two configurations is possible 
via mixing of the positive- and negative-parity levels (right-hand side). Octupolc 
deformation (3% of the mean field rises in the crossing region. In this figure, the 
crossing between the two levels with the asymptotic quantum numbers [321 1] and 
[200tj] is illustrated as an example. The two levels satisfy the selection rule, An^ = 2 
and A A = 1, for the matrix elements of the non-axial octupole operator r^Y^i, so 
that the mixing between them takes place mainly through the r 3 Y3i component of 
the mean field. 
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Figure 6: Plot of the excitation energies versus angular-momentum for the SD ro- 
tational bands in 32 S, 36 S, and 50 S calculated by means of the cranked SHF method. 
Results obtained with the use of the SIII, SkM*, and SLy4 interactions are plotted by 
solid, dashed, and dotted curves, respectively. Their shape evolutions as functions 
of angular momentum / in the (/?2>7) plane are displayed in the upper portions. 
The (52 values decrease with increasing /. 
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Figure 7: The angular momenta / and the moments of inertia J are plotted as 
functions of rotational frequency u; ro t for the SD rotational bands in 32 S, 36 S, and 
50 S. The SLy4 interaction is used. Values of the kinematical and dynamical moments 
of inertia, = I/uj TOt and = dI/duj ro t, are plotted in unit of h 2 /MeV by 
solid and dashed curves, respectively. For reference, the rigid-body moments of 
inertia J7ri g = m J p(r)(y 2 + z 2 )dr evaluated with the calculated density p{r) are 
also indicated. 
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Figure 8: Single-particle energy diagrams (Routhians) for the SD bands in 32 S, 36 S, 
and 50 S, plotted as functions of rotational frequency uj ro t- The left(right)-hand side 
displays those for protons(neutrons). The levels associated with the gg/ 2 and f 7 / 2 
shells are drawn by thick-solid and thick-dashed lines, respectively. Other occupied 
levels associated with the sd and fp shells are drawn by thin-solid and thin-dashed 
lines, respectively. Unoccupied levels are drawn by thin-dotted lines. Numbers 
indicate the Fermi surfaces and total numbers of single-particle states below them. 
The result calculated with SLy4 is shown here, but we obtained similar results also 
with SIII and SkM*. 
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Figure 9: Deformation energy curves (measured from energies at = 0) as func- 
tions of the octupole deformation parameters I3^ m {m = 0, 1, 2, 3), calculated for 32 S, 
36 S, and 50 S, by means of the constrained HF procedure with the use of the SIII, 
SkM* and SLy4 interactions. The quadrupole deformation parameters are fixed at 
the equilibrium value of in each nucleus and 7 = 0. One of the (5^ m (m = 0, 1, 2, 3) 
is varied while the other /33 m 's are fixed to zero. The deformation parameters f3% 
and are defined in terms of the expectation values of the octupole operators 
(see Ref. |18| for their explicit expressions). 
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Figure 10: Left-hand side: Density distributions in the (y,x)- and (z,x)- planes 
of the SD band at / = in 32 S, 36 S, and 50 S, calculated with the use of the SLy4 
interaction. Neutron (proton) equi-density lines with 50% and 1% of the central 
density are shown by dashed (solid) lines (the inner and outer lines correspond to 
the 50% and 1% lines, respectively). Right-hand side: Same as the left-hand side 
but for I = 20, 22, 28 for 32 S, 36 S, and 50 S, respectively. 
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